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Abstract. Several problems at the interface between the field-theoretical description of the 
Casimir effect and experiments on measuring the Casimir force are discussed. One of these 
problems is connected with the definition of the Casimir free energy in ideal metal rectangular 
boxes satisfying the general physical requirements. It is shown that the consideration of 
rectangular boxes with a partition (piston) does not negate the previously known results 
obtained for boxes without a piston. Both sets of results are found to be in mutual agreement. 
Another problem is related to the use of the proximity force approximation for the interpretation 
of the experimental data and to the search of analytical results beyond the PFA based on the 
first principles of quantum field theory. Next, we discuss concepts of experimental precision 
and of the measure of agreement between experiment and theory. The fundamental difference 
between these two concepts is clarified. Finally, recent approach to the thermal Casimir force 
taking screening effects into account is applied to real metals. It is shown that this approach is 
thermodynamically and experimentally inconsistent. The physical reasons of this inconsistency 
are connected with the violation of thermal equilibrium which is the basic applicability condition 
of the Lifshitz theory. 



1. Introduction 

For almost 50 years since Casimir's discovery [1], theory existed independent of rare experiments 
[2,3]. A lot of theoretical work has been done during this period. However, the relationship 
to reality of some theoretical models, such as an ideal metal spherical shell, an ideal metal 
rectangular box, a dielectric ball etc., remains unclear up to now. In the last ten years 
scientific investigations in the field of the Casimir effect have experienced an interaction between 
experiment and theory. This has revealed that the application of some basic theories to real 
experimental situations is highly nontrivial and even leads to communication difficulties between 
theorists and experimentallists. 

In this paper we summarize some experience of the interaction between "high theory" and 
real world experimental details in the last ten years. Different points of view are considered 
on such problems as agreement between experiment and theory and applicability of some ideal 
models and approximate methods in real experimental situations. It is shown that in some cases 
confusion arises from an inadequate use of terminology. 

In Sec. 2 we discuss an old problem of the thermal Casimir force in ideal metal rectangular 
boxes and suggest a new solution which satisfies general physical criteria. It is shown that the 
case of an isolated box is independent of a box with a partition (piston). The results for the 
Casimir force obtained for each of these configurations are in mutual agreement. 
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Section 3 briefly reviews the proximity force approximation including its justification from 
the first principles of quantum field theory and experimental applications. In this respect a new 
representation for the Casimir energy in terms of the functional determinants and scattering 
matrices is considered. The application of this representation to real material bodies is still 
problematic. 

In Sec. 4 we consider the problem of the reliability of experiments. It is underlined that the 
experimental error is an independent characteristic of the experimental precision which should 
not be confused with the measure of agreement between experiment and theory. 

Section 5 is devoted to the comparison between experiment and theory in the measurements 
of the Casimir force. In this respect different approaches to the theoretical description of the 
Casimir force between real metals are compared with the most precise indirect measurement of 
the Casimir pressure between two parallel plates by means of micromechanical torsional oscillator 
[4,5]. Special attention is paid to uncertainties which might be introduced in the computations 
due to deviations of the tabulated optical data from the data particular to the metallic films 
actually used. 

In Sec. 6 we consider a recent theoretical approach to the thermal Casimir force taking into 
account the screening effects and diffusion currents. We apply this approach to the case of 
real metals and analyze its consistency with the principles of thermodynamics and experimental 
data. Specifically we show that for metals with perfect crystal lattices the inclusion of screening 
effects results in violation of the Nernst heat theorem. The experimental data of the experiment 
[4,5] exclude this approach at a 99.9% confidence level. 

Section 7 contains our conclusions and discussion. 



2. Thermal Casimir force in ideal metal rectangular boxes 

Ideal metal rectangular boxes were first considered by Lukosz [6], Mamayev and Trunov [7,8] and 
Ambj0rm and Wolfram [9]. This configuration attracted much attention because it was found 
that the electromagnetic Casimir force in rectangular boxes can be both attractive and repulsive 
depending on the ratio of sides a x , a y and a z along the x, y and z axes. The nonrenormalized 
Casimir energy of the box is equal to (for simplicity we consider the massless scalar field with 
Dirichlet boundary conditions) 
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The regularization of ([T]) can be performed, e.g., using the Epstein zeta function or the cut-off 
method [10]. The latter permits to find the geometric structure of infinities contained in ([T]). 

To do so, one replaces Eq(ci x , a y , a z ) from ([I]) with EQ S \a x ,a y ,a z ) by introducing the cut-off 
function 

f(5uo nlp ) = e~ s ^ (3) 

under the sign of summation in (pQ). After the repeated application of the Abel-Plana formula 

[10] to Eq^ (a x ,a y , a z ) one finds that there are three different types of divergent quantities in the 
limit 8 — > 0, Ii, I2 and I3 of order 5~ 4 , 5 -3 and 5~ 2 , respectively. Then, the finite, renormalized, 
Casimir energy can be defined as 
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Here, Ik (k = 1, 2, 3) are the counter terms having the following geometrical structure: 
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A similar situation takes place for the electromagnetic field, where the renormalized Casimir 
energy, E T ^ m) also takes the form of (with Eq replaced for Eq 111 ) and 

If m = 2/i, I| m = 0, I% m = -2/3. (6) 

It is seen that in both cases the counter terms are proportional to the volume of the box 
V to the area of box surface and to the sum of sides. 

In the last few years the configuration of a rectangular box with so-called movable partition 
(piston) has attracted much attention [11—14]. This means that the piston can have any fixed 
position parallel to the two opposite faces of the box (the configuration where the piston is not 
fixed and may slide between the opposite faces is in fact a nonequilibrium case). Let the piston 
be parallel to the plane xy and have an equation z = a z \ < a z . In this case our box is divided 
into the two boxes a x x a y x a z \ and a x x a y x (a z — a z \). Calculating the sum of the regularized 
Casimir energies 

Eq (a x , a y , a z \) + E (a x , fly, CL Z Cl Z l 

), (7) 

one finds that the contribution from the singular terms of the form of © does not depend on 
the position of the piston a z \. This leads to a finite force acting on the piston 
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This force is well defined and does not require the renormalization procedure (jlj). 

In both scalar and electromagnetic cases the force acting on the piston attracts it to the 
nearest face of the box. On this ground the existence of the Casimir repulsion in cubes in the 
electromagnetic case was cosidered doubtful [12]. Specifically, it was claimed [12,13] that the 
definition of the pressure acting on a cube face requires elastic deformations of single bodies 
treated as perfect. The attraction (or repulsion for a piston with Neumann boundary conditions 
[15]) of a piston to the nearest face of the box does not, however, negate the Casimir repulsion 
for boxes without a piston that have some appropriate ratio of a x , a y and a z . The point is that 
the cases with an empty space outside the box and that with another section of the larger box 
outside the piston are physically quite different. In the first case the vacuum energy outside 
the box does not depend on a x , a y and a z and there is no force acting on the box from the 
outside. Whereas in the second case there is an extra section of the larger box outside the 
piston which gives rise to the additional force acting on it. In fact one need not admit elastic 
deformations to define a force and a pressure in static configurations. This is simply done using 
the principle of virtual work and virtual displacements through real forces [16,17]. In addition, 
from a thermodynamic point of view any equilibrium system can be characterized by the free 
energy (energy if the temperature is equal to zero) and the respective pressure [18] 
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From this point of view it would be illogical to admit consideration of the force acting on a 
piston, but exclude from consideration forces acting on the faces of a box where this piston 
serves as a partition. 

In this respect it seems important to provide a finite definition of the Casimir free energy 
in ideal metal rectangular boxes satisfying general physical requirements. The first calculations 



on this subject [9] resulted in a divergent free energy after removing the regularization. More 
recent results appear to be either infinite [19] or ambiguous [20]. Paper [21] reconsidered the 
derivation of the Casimir free energy in rectangular boxes using zeta functional regularization. 
However, the used formalism does not include all necessary subtractions. 

The following definition of the Casimir free energy in rectangular boxes suggests itself 
[12,13,21] 
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This expression is finite. However, it cannot be considered as physically satisfactory. The 
problem is that at high temperature the thermal correction ()10p behaves as [22] 
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where a\ = — Vn 2 /90, 02,3 = tt2,3(a x , a y , a z ) can be expressed in terms of the heat kernel 
coefficients and 04 = const. Then at high temperature AtJ~o contains terms of quantum 
origin which increase with the increase of temperature. In the general case, these terms lead to 
respective forces acting on the box faces which increase with the increase of the sides of the box. 
Such paradoxical properties are physically unacceptable. Because of this it was suggested [23] 
to define the physical Casimir free energy of the box as 
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With this definition, the respective Casimir forces acting on the box faces go to zero when all 
the box sides a x , a y , a z go to infinity in agreement with physical intuition. 

The physical meaning of all three subtractions made on the right-hand side of (|12p can be 
clearly understood. The first term is actually the contribution of the blackbody radiation in 
the volume of the box. This is seen from the fact that the free energy density of the blackbody 
radiation in empty space is given by 

f d 3 k , / _Mhi\ ir 2 (k B T) 4 , . 

where for the electromagnetic case = 2/^. 

For the scalar Casimir effect in a rectangular box with sides a x x a y x a z the asymptotic 
behavior of AxJ~o at high T was investigated in [23] with the result 
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In the electromagnetic case the following values of these coefficients were obtained: 



7T 



al m = 0, af* = —(a x + a y + a z ). (15) 

This demonstrates that the geometric structures of all three terms subtracted in (|12|) are precisely 
the same as the terms subtracted in (jl]) to make the Casimir energy finite at zero temperature. 
Because of this, the subtraction procedure in (|12j) can be interpreted as the additional (finite) 
renormalization of the same geometric parameters as were renormalized at zero temperature to 
make the Casimir energy of the box finite. 



The simplest application of the final expression for the physical Casimir free energy (I12D is 
the case of two plane parallel plates. It is easily seen that in this configuration «2 = «3 = and 
one is left with only a subtraction of the free energy of the blackbody radiation in the volume 
between the plates V = aS, where S is the infinite plate area. This leads to the well known 
result [10,24] for the electromagnetic Casimir free energy per unit area of the plates 
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where t = T e g/T, and the effective temperature is defined as ksT^ = hc/(2a). In particular, at 
T <C T c ff one obtains 
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where the last contribution on the right-hand side originates from the subtraction of the 
blackbody radiation. We emphasize that only this term contributes to the thermal correction 
to the electromagnetic Casimir pressure at low temperatures (short separations) 
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Equation (|12p solves the long-standing problem on the calculation of the physical Casimir 
free energies and pressures in rectangular boxes of any size. A few examples for both the scalar 
and electromagnetic Casimir effect are considered in [23]. Here we present the computational 
results for the electromagnetic free energy in a cube and for the respective Casimir force 
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acting on the opposite cube faces. 

In Fig. 1(a) we plot the electromagnetic Casimir free energy in a cube as a function of a at 
T = 300 K (solid line). In the same figure the Casimir energy at T = is shown by the dashed 
line. As is seen in this figure, the electromagnetic Casimir free energy decreases with the increase 
of separation. At large separations T approaches a constant. In Fig. 1(b) the electromagnetic 
Casimir free energy is shown as a function of temperature for a cube with a = 2 fim. The free 




Figure 1. The electromagnetic Casimir free energy for a cube as a function of (a) size a at 
T = 300 K (solid line; the dashed line shows the energy at T = 0) and (b) temperature at 
a = 2 fim. 
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Figure 2. The electromagnetic Casimir force between the opposite faces of a cube as a function 
of (a) size a at T = 300 K (solid line; the dashed line shows the force at T = 0) and (b) 
temperature at a = 2 fim. 



energy decreases with the increase of T. At high temperatures T demonstrates the classical 
limit. The respective thermal electromagnetic Casimir force at T = 300 K, as a function of 
a, is shown in Fig. 2(a) by the solid line. It is positive (i.e., repulsive) for cubes of any size. 
Thus, thermal effects for cubes in the electromagnetic case increase the strength of the Casimir 
repulsion. The dashed line in Fig. 2(a) shows the electromagnetic Casimir force at T = as a 
function of a. This force is given by 

- (20) 

i.e., it is always repulsive. Fig. 2(b) demonstrates the electromagnetic Casimir force in a cube 
of size a = 2 fim as a function of temperature. It is seen that the force increases with increasing 
temperature. 

Note that the results presented differ from those found in [21] where the terms of order 
(/cgT) 4 and of lower orders in the Casimir free energy were obtained in the high-temperature 
regime. This is explained by the fact that the authors of [21] did not make subtractions of the 
contributions from the blackbody radiation and of the terms proportional to the box surface 
area and to the sum of its sides. 

The thermal correction to the Casimir energy and force acting on a piston were investigated 
in [13] for the scalar field with Dirichlet or Neumann boundary conditions using the definition 
(|10p . The electromagnetic Casimir free energy and force acting on a piston were found in the 
case of ideal metal rectangular boxes and cavities with the general cross section [13]. In the limit 
of low temperatures the thermal correction to the Casimir force on a piston was shown to be 
exponentially small. In the case of medium temperature a x <C hc/ftsT) <C a y ,a z the authors 
of [13] obtained terms of order (A^T) 4 and of order (kBT) 2 in the electromagnetic Casimir free 
energy. In the scalar Casimir free energy, a term of order (fc^T) 3 was also obtained. This 
results in the contribution to the force acting on a piston which increases with the increase of 
the temperature, depends on h and c and does not depend on the position of the piston. The 
scalar and electromagnetic thermal Casimir forces acting on a piston were also considered on 
the basis of equation ([TO]) in [25]. 

The same results for the thermal correction to the Casimir force acting on a piston are 
obtained if the free energy is defined in accordance with equation (|12|) . This is because the 
contribution of blackbody radiation to the energy of the entire box is equal to 



— a x a y a z \ fob — a x a y (a z — a z \)fob = —a x a y a z fob, (21) 
i.e., it does not depend on the position of the piston. This is also true for terms of order (fc^T) 3 



and {ksT) 2 which are proportional to the surface area of each section of the box and to the sum 
of its sides. 

The above results were obtained for rectangular boxes with the Dirichlet boundary conditions 
(scalar case) and for ideal metal boxes (electromagnetic case). In the same way, as for 
zero temperature, the consideration of the thermal Casimir effect in rectangular boxes has to 
incorporate real material properties of the boundary surfaces. Till now this problem has not 
been conclusively solved. 



3. Functional determinants and the justification of the proximity force 
approximation 

The proximity force approximation [26] provides an important bridge between experiment and 
theory. Experimentally it is hard to use the configuration of two parallel plates. Because of this, 
most of experiments use the configuration of a sphere above a plate for which, even in the ideal 
metal case, the exact results for the electromagnetic Casimir force are not available. According 
to the proximity force approximation (PFA), the interaction energy between two curved surfaces 
Si and T,2 can be approximately calculated by replacing the small curved surface elements with 
respective plane plates. If the interaction energy between the opposite plane parallel elements 
is notated as E(z) (where z is the separation distance), the interaction energy and force are 
approximately represented as 

U(a) = jf E(z)da, F(a) = (22) 

For the configuration of an ideal metal sphere of radius R at a separation a above an ideal metal 
plane ([22]) results in 

-^pfa ( a ) = 2*RE{a) = -^J. (23) 
For an ideal metal cylinder above an ideal metal plate the PFA leads to 
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Equations (|22p - (j24jl are the approximate ones. They are applicable only at short separations 
between the surfaces. Thus, ([23]) and ([2*31) work well only at a -C R. 

In many papers the PFA (I22p is applied in a region where it is not applicable, for example 
at a = R/2. The obtained large deviations of the PFA result from the exact result are then 
considered as a "violation of the PFA". Such formulations are in fact misleading. The PFA 
gives only the main contribution to the force under some conditions. Specifically, it would be 
meaningless to calculate the integral in (|22|) up to higher orders in the related small parameter 
with the aim of obtaining a more exact result. What is really meaningful is the search of an 
exact analytical representation for the Casimir force in configurations where only the PFA result 
is so far available. 

In the last few years the finite representation for the Casimir energy for two separated bodies 
A and B in terms of the functional determinants was obtained. In this representation the 
Casimir energy can be written in the form [27,28] 
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for a body A and B, respectively. The latter is widely used in light scattering theory, where it is 
the basic object for expressing the properties of the scatterers [29]. Using such a representation, 
in [28] the analytic results for the electromagnetic Casimir energy for an ideal metal cylinder 
above an ideal metal plane were obtained. Eventually, the result is expressed through the 
determinant of an infinite matrix with elements given in terms of the Bessel functions. The 
analytic asymptotic behavior of the exact Casimir energy at short separations was found in [30] . 
It results in the following expression for the Casimir force at a <C R: 



The PFA result (j24j) in this case matches with the first term on the right-hand side of (|26|) , 

Equation (|26|) is very important. It demonstrates that the relative error of the electromagnetic 
Casimir force between a cylinder and a plate calculated using the PFA is equal to 0.2886 a/R. 
Thus, for typical experimental parameters of R = 100 /im and a = 100 nm this error is 
approximately equal to only 0.03%. 

For a sphere above a plate made of ideal metals the exact analytic solution in the 
electromagnetic case has not yet been obtained. The scalar Casimir energy for a sphere above a 
plate was found in [30,31]. The scalar Casimir energies for both a sphere and a cylinder above 
a plate have also been computed numerically using the wordline algorithms [32,33], but it was 
noted that the Casimir energies for the Dirichlet scalar field should not be taken as an estimate 
for those in the electromagnetic case. For an ideal metal sphere above an ideal metal plane a 
correction of order a/R beyond the PFA was computed numerically in [34] for a/R> 0.075 and 
in [35] for a/R > 0.15. In both cases the extrapolation of the obtained results to smaller a/R 
leads to a coefficient near a/R approximately equal to 1.4. 

In addition, the validity of the PFA for a sphere above a plate has been estimated 
experimentally [36] and the error introduced from the use of this approximation was shown 
to be less than a/R. This is in disagreement with the extrapolations made in [34,35]. To solve 
this contradiction, it is desirable to find the analytical form of the first correction beyond the 
PFA for a sphere above a plane, like in (|26p for the cylinder-plane configuration. 

In fact the representation (|25[) provides a far-reaching generalization of the Lifshitz formula. 
From conceptual point of view it can be applied not only to ideal metals, but to real materials 
as well. The problem, however, is to find the matrix elements of the T-matrix operator which 
would take proper account of both geometric shape and material properties of the test bodies 
used in the experimental situation. 

4. The experimental error and reliability of experiments 

The concept of the experimental error is often confused with the theoretical error and with 
the measure of agreement between experiment and theory. However, when we deal with an 
independent measurement, the experimental error has nothing to do with any theory of the 
measured quantity. The independent measurement of the Casimir force or its gradient does not 
use any theory of the Casimir effect. Thus, the experiments [4,5,37-42] are independent in this 
respect. In other experiments (in [43], for instance) the measurement data are fitted to some 
theoretical expression for the Casimir force. Such kind of measurements are not independent 
and we do not consider them below. 

Some papers arrive to theoretical conclusions which are inconsistent with the measurement 
data. This is sometimes surrounded by the statement that the measurements might be not as 
precise as indicated by the authors. It is our opinion that such statements made without an 
indication of any specific cause are inappropriate. Both random, A rand i Jiexpt (a), and systematic, 
A syst i ?expt (a), experimental errors in the Casimir force measurements are found using the 
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Figure 3. The total absolute experimental error of the Casimir pressure measurements [4,5] (the 
solid line), the random error (the long-dashed line), and the systematic error (the short-dashed 
line) are shown as functions of separation. 



rigorous statistical procedures. They can be combined to find the total experimental error 

A tot F cxpt (a) = q p {r) [A rand F cxpt (a) + A syst F cxpt (a)l . (27) 

Here, qp{r) determined at j3 = 0.95 (i.e., at 95% confidence level) varies between 0.71 and 0.81 
depending of the value of the quantity r = A sy F exp (a)/s^(a), where sp(a) is the variance of 
the mean measured quantity [44]. In fact there is no arbitrariness in the determination of the 
total experimental error which is the ultimate characteristic of the precision of the measurements. 
The most valuable esperiments are marked by a negligible role of the random error. For such 
experiments 

A tot F expt (a) « A syst F cxpt (a). (28) 

For today there is only one indirect measurement of the Casimir pressure between Au coated 
plates by means of micromechanical torsional oscillator satisfying this condition [4,5]. The 
total experimental error in this measurement at shortest separations is as small as 0.2% of 
the measured Casimir pressure. We stress once again that this error is unrelated to much 
larger errors inherent to theoretical computations on the basis of the Lifshitz theory or to the 
measure of agreement between experiment and theory. This is just the resulting error with 
which the experimental data are taken. As an example, the total absolute experimental error 
in the experiment on measuring the Casimir pressure by means of a micromechanical torsional 
oscillator [4,5] is shown in Fig. 3 as a function of separation (the solid line). The long-dashed and 
short-dashed lines show the random and systematic errors, respectively. As a result, the relative 
total experimental error 5 tot P cxpt (a) = A tot P expt (a)/|P expt (a)| varies from 0.19% at a = 162 nm 
to 0.9% at a = 400 nm, and to 9.0% at a = 746 nm. 

Sometimes the experimental precision can be questioned if there are some doubts in the 
calibration procedures used. For example, the electrostatic calibration is of prime importance in 
the independent measurements of the Casimir force. Specifically, it is usually carefully verified 
that the residual potential between the grounded test bodies does not depend on separation 
where the measurements of the electric force are performed. Recently it was claimed that 
the residual potential Vq from the electrostatic calibration in the sphere-plate configuration is 
separation dependent [45]. The authors used an Au-coated sphere of 30.9 mm radius above an 
Au coated plate. On the basis of these measurements a reanalysis of the independence of Vq 
on separation in the earlier measurements of the Casimir force by means of an atomic force 
microscope and a micromachined oscillator was invited. The results [45] are, however, not 
directly relevant to the earlier measurements. The point is that the radius of the sphere used in 
[45] is a factor of 300 larger than in the earlier precision measurements of the Casimir force. It 
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Figure 4. The crosses show the measured mean Casimir pressures together with the absolute 
errors in the separation and pressure as a function of the separation, (a) The theoretical Casimir 
pressures computed using the generalized plasma-like model and the optical data extrapolated 
by the Drude model are shown by the light-gray and dark-gray bands, respectively, (b) The 
theoretical Casimir pressures computed using different sets of optical data available in the 
literature versus separation are shown as the dark-gray band. 



is well known that for large test bodies (i.e., large interaction areas) there are large variations 
of electric forces due to deviations of the mechanically polished and ground lens surface from 
perfect spherical shape [46]. 

5. Comparison between experiment and theory 

Experiment is the supreme arbiter in physics. Because of this, the comparison between 
experiment and theory is a painful point for those theories that are found to be experimentally 
inconsistent. It happens that in such cases both the experimental data and the methods of 
comparison are questioned. The Casimir force is a strongly nonlinear function of the separation 
distance. As a consequence, such global characteristics of the agreement between experiment 
and theory as the root-mean-square deviation were found to be inadequate [47]. In the last 
few years two local methods on how to compare experiment with theory in the Casimir force 
measurements were elaborated and successfully applied. Within the first method [4,48,49], the 
experimental data are represented as crosses with arms determined by the total experimental 
errors in the measurement of separation and a related quantity (the force, the pressure or the 
frequency shift) determined at some chosen confidence level. In the same figure, one should plot 
the theoretical band whose width is equal to the total theoretical error determined at the same 
confidence as the experimental errors. The overlap (or its absence) of the experimental crosses 
and the theoretical band can be used to make a conclusion on the consistency or inconsistency 
between experiment and theory. 

In Fig. 4 the first method of comparison between experiment and theory is illustrated on the 
measurement data by Decca et al. [4,5] discussed in Sec. 3. The light-gray band in Fig. 4(a) 
shows the theoretical results computed using the Lifshitz theory combined with the generalized 
plasma-like dielectric permittivity [50,51] 
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Here, uj p is the plasma frequency, Uj ^ are the frequencies of the oscillators describing core 
electrons, fj are the oscillator strengths and jj are the relaxation parameters. The dark-gray 
band in Fig. 4(a) is computed by the same Lifshitz theory using the tabulated optical data for 
Au [52] extrapolated to low frequencies by means of the Drude model [53-55] 
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Here c(i£) is the conductivity. It is connected with the dc conductivity by the equation 
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Note that the plasma frequency and the dc conductivity are expressed as [56] 
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where e and m are the charge and the mass of an electron, n is the charge carrier density and /x is 
their mobility. As is seen in Fig. 4(a), the experimental data shown as crosses (the experimental 
errors are determined at a 95% confidence level) are consistent with the theoretical approach 
using the generalized plasma-like permittivity. The Drude model approach is excluded at a 95% 
confidence level. In Fig. 4(b) the same experimental data are reproduced and compared with the 
Drude model approach using all sets of optical data available in the literature [57]. As is seen 
in Fig. 4(b), the use of optical data alternative to [52] makes the disagreement deeper between 
the experimental data and the Drude model approach. In Fig. 4(a,b) the comparison between 
experiment and theory is performed within the separation region from 500 to 600 nm. However, 
exactly the same conclusions follow over the entire measurement range in this experiment from 
160 to 750 nm. 

In the second method for the comparison of experiment and theory in the Casimir 
force measurements [40,48,58], the differences between the theoretical and mean experimental 
quantity, for instance, P theor (a) — P expt (a), are plotted as dots. In the same figure the borders of 
the confidence intervals [— Sp(a), Hp(a)] for this difference at a chosen confidence level (usually 
95%) are plotted as the function of separation. If no less than 95% of the dots representing 
the above differences belong to the confidence interval the theoretical approach is consistent 
with the data. Alternatively, if almost all the dots are outside the confidence interval, the 
theoretical approach is excluded by the data at a 95% confidence level. In Fig. 5 we illustrate 
the second method for the comparison of experiemnt with theory using the experimental data 
of the same measurements [4,5]. In Fig. 5(a) the theoretical approach using the generalized 
plasma-like permittivity (|29|) is compared with the data. It is seen that all dots are inside the 
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Figure 5. The differences of the theoretical and the mean experimental Casimir pressures 
between the two Au plates versus separation are shown as dots. The theoretical results are 
calculated using the Lifshitz theory at room temparature using (a) the generalized plasma-like 
model and (b) the Drude model approach. The solid lines indicate the boundaries of the 95% 
confidence intervals. The dashed line indicates the boundary of the 99.9% confidence intervals. 



error bars. Thus, this approach is consistent with the data. In Fig. 5(b) the same data are 
compared with the theoretical approach using the tabulated optical data extrapolated by the 
Drude model. The solid and dashed line represent the borders of 95% and 99.9% confidence 
intervals, respectively. As is seen in Fig. 5(b), the Drude model approach is experimentally 
excluded at a 95% confidence level within the entire measurement range from 160 to 750 nm. 
Within a more narrow measurement range from 210 to 620 nm the Drude model approach is 
excluded at a 99.9% confidence level. 

If the theoretical approach is experimentally consistent [see Fig. 5(a)] the quantity Hp/|P expt |, 
determined at a 95% confidence level, can be used as the quantitative measure of agreement 
between experiment and theory. Thus, at a = 162 nm this measure is equal to 1.9%. It decreases 
to 1.4% at a = 300 nm and than gradually increases up to 9.7% at a = 745 nm. It is evident that 
at the shortest separation the agreement between experiment and theory is almost an order of 
magnitude worse than the total experimental error equal to only 0.19%. This is explained by 
large theoretical errors which dominate in the determination of Hp at the shortest separations. 

The above explanations aim to make absolutely clear that the calculation of errors and the 
comparison between experiment and theory is not an arbitrary, but a rigorously determined 
procedure. Recently, the measurement data of the experiment [4,5] was independently 
reanalyzed in [59] with the conclusion: "The data rule out the Drude approach. . . , while they 
are consistent with the plasma-model approach. . ." 



6. Attempt to account for screening effects 

Recently, the above discussed problems of the Drude model approach in application to 
real metals, and related problems arising for dielectric and semiconductor materials [60-64], 
motivated an attempt to modify the reflection coefficients in the Lifshitz formula by including 
the screening effects and diffusion currents [65,66]. The modified reflection coefficients for the 
transverse magnetic and transverse electric modes were obtained through use of Boltzmann 
transport equation which takes into account not only the standard drift current j, but also the 
diffusion current eDVn, where D is the diffusion coefficient and Vn is the gradient of the charge 
carrier density [66]. The transverse magnetic coefficient takes the form 



?"TM(i£,fcj 



my. 



(33) 



where k± is the projection of the wave vector in the plane of the plates, u = i£ is the imaginary 
frequency and the following notations are introduced 



2 1 2 i £ 

cr 



e(i£) m) - e(i£) 



1/2 



(34) 



In this equation, 1/k is the screening length and the dielectric permittivity of core electrons 
is defined in (|29p . The transverse electric coefficient is given by the standard expression 



k 



q + k' 



(35) 



as is used in the Drude model approach. 



The paper [66] claims the application of the above approach to intrinsic semiconductors only. 
It uses a specific Debye-Hiickel expression for the screening length 




i = ^ =Rm= Jm^. ( 36) 

K kdh V 47re z n 

This expression is applicable to particles obeying the Maxwell-Boltzmann statistics. It is 
obtained from the general representation for the screening length [67] 

(37) 

if one uses the expression (I32D for the dc conductivity and Einstein's relation [56,67] 

^ = ^ («) 
fj, \e\ 

valid in the case of Maxwell-Boltzmann statistics. In the limiting case £ — > the reflection 
coefficient ([3*3*]) coincides with that obtained in [65] . 

However, the application region of the reflection coefficients (|33|) . (]34|) with the Debye-Hiickel 
screening length (|36p cannot be restricted to only intrinsic semiconductors. These coefficients 
should be applicable to all materials where the density of charge carriers is not too large so 
that they are described by Maxwell-Boltzmann statistics. This means that in the framework 
of the proposed approach it is legal to apply ([33]) - ([36]) to doped semiconductors with dopant 
concentration below critical and to solids with ionic conductivity etc. 

Here, we consider the application of this approach to metallic plates. Metals and 
semiconductors of metallic type are characterized by rather high concentration of charge carriers 
which obey the quantum Fermi-Dirac statistics. The general transport equation, however, is 
equally applicable to classical and quantum systems. The only difference one should take into 
account is the type of statistics. Substituting Einstein's relation, valid in the case of Fermi-Dirac 
statistics [56,67] 

- = ff- ™ 

[i 3|e| 

where Ep = Tiu) p is the Fermi energy, into (j37[) . one arrives to the following expression for the 
Thomas-Fermi screening length [67] 
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With this definition of the parameter k, it is legal to apply equations ([33]) - ([35] ) to metals. 

Now we consider two thick metallic plates separated by a distance a at temperature T in 
thermal equilibrium. Under these conditions the Casimir free energy per unit area of the plates 
is given by the Lifshitz formula [68]. Let us assume that the reflection coefficients (I33p - (l35p . 
(|4Up can be substituted into this formula. Then in terms of dimensionless variables y = 2aq, 
£ = i/uJ c = 2a£ / c one obtains 

7 rr~\ OO n 

?(a,T) = -^J2 / ydy{ln\l-4 M (i( h y)e-y] + In [l - r^(i( h y) e""] } , (41) 



where Q = AirakBTl/{Tic) are the dimensionless Matsubara frequencies and a prime near the 
summation sign adds a multiple 1/2 to the term with 1 = 0. In terms of the dimensionless 
variables the reflection coefficient (133)) takes the form 



where 
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Note that all dielectric permittivities here are functions of iui c C- Below we do not use the explicit 
expression for the reflection coefficient f^E^dy) because it coincides with the standard one, as 
defined in the Drude model approach, and considered in detail in [69]. 

Let us determine the behavior of the Casimir free energy (I4ip at low temperature. For all 
metals the screening length (I40p is very small. As a result, at any reasonable separation distance 
between the plates, the dimensionless parameter n a defined in (|43f) is very large and the inverse 
quantity can be used as a small parameter 



2oktf = K a > 1, 
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(44) 



Expanding the reflection coefficient (|42l) up to the first power of the parameter (3 a one obtains 

r TM (iC, V) = r T M(iC, y) ~W a Z + 0(J5 2 a ), (45) 
je{e -e) 3 y(y 2 
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where rTM(iCi2/) is the standard TM reflection coefficient calculated with the dielectric 
permittivity e(iu; c C) [it is given by (|42p with the third term in both numerator and denominator 
omitted]. From (|45p one arrives at 
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Now we substitute (|46|) and the respective known expression for the TE contribution [69] 
into (|4ip . Calculating the sum with the help of the Abel-Plana formula, we obtain in perfect 
analogy to [69] 
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where J-^'(a,T) is determined by equation (17) in [69]. It goes to zero together with its 
derivative with respect to temperature when T — > 0. The quantity J 7 ^ (a, T) originates from 
the second contribution on the right-hand side of (I45p . It is easily seen that j^(a,T) = 
E^\a) + 0(T 3 /T e 3 ff ) at low T. The Casimir free energy J r gp (a, T) is defined by substituting the 
dielectric permittivity (I29p of the generalized plasma-like model into the Lifshitz formula. It 
was found in [50,51] and the respective thermal correction was shown to be of order (T/T e fj) 3 
when T — > 0. The TE reflection coefficient at zero frequency entering (|47f) is given by 
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As a result, calculating the Casimir entropy 
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with the use of (|47p and considering the limiting case of zero temperature, one arrives at 
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in violation of the Nernst heat theorem. This result is obtained for metals with perfect crystal 
lattices. In the presence of impurities the Casimir entropy abruptly jumps to zero at T < 10" 3 K 

Thus, the modified reflection coefficients taking the screening effects into account lead to a 
violation of the Nernst heat theorem for metals with perfect crystal lattices in the same way 
as the standard Drude model approach. Because of this, the theoretical approach using such 
reflection coefficients is thermodynamically inconsistent. 

Now we briefly compare the theoretical predictions, following from the use of reflection 
coefficients rpM and ttEj with the measurement data of the most precise experiment by means 
of micromachined torsional oscillator [4,5]. This experiment was already discussed in Sec. 5. In 
Fig. 6(a) the experimental data for the Casimir pressure between two Au plates are shown as 
crosses with the absolute errors determined at a 95% confidence level. The solid line presents the 
computational results for P(a, T) = —d!F(a : T) /da using the Lifshitz formula and the generalized 
plasma- like dielectric permittivity (|29p . The parameters of oscillators for Au were determined 
in [5] with high precision. The dashed line was computed using the Lifshitz formula for 
P mod (a,T) = —dJ-(a,T)j da with the reflection coefficients ttm,te taking the screening effects 
into account. As is seen in Fig. 6(a), the theoretical approach taking into account the Thomas- 
Fermi screening length is experimentally excluded at a 95% confidence level over the separation 
region from 500 to 600 nm. The same conclusion follows within the entire measurement range 
from 160 to 750 nm. 
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Figure 6. (a) The crosses show the measured mean Casimir pressures together with the absolute 
errors as a function of the separation. The theoretical Casimir pressures computed using the 
generalized plasma-like model and the approach including the screening effects are shown as solid 
and dashed lines, respectively, (b) Differences of the theoretical Casimir pressures computed with 
inclusion of the screening effects and the mean experimental Casimir pressures versus separation 
are shown as dots. The 95% and 99.9% confidence intervals are shown as the solid and dashed 
lines, respectively. 



Fig. 6(a) illustrates the first method for the comparison between experiment and theory 
in Casimir force measurements discussed in Sec. 5. In Fig. 6(b) the second method for the 
comparison of experiment and theory is illustrated. Here, the differences between the theoretical 
Casimir pressures computed with inclusion of the screening effects and the mean experimental 
pressures are shown as dots. The solid line indicates the borders of 95% confidence intervals. 
Dots are outside the confidence interval [Sp(a),'Ep(a)] over the entire measurement range 
from 160 to 750 nm. In the same figure, the dashed line shows the borders of 99.9% confidence 
intervals. As is seen in Fig. 6(b), dots are outside of this confidence interval within the separation 
region from 160 to 640 nm. Thus, within this region of separations the theoretical approach 
taking the screening effects into account [66] is experimentally excluded at a 99.9% confidence 
level. 

The physical reasons why the inclusion of the screening effects into the Lifshitz theory is 
thermodynamically and experimentally inconsistent can be understood as follows. The Lifshitz 
theory is formulated for systems in thermal equilibrium. As was indicated in [71], the drift 
current of conduction electrons leads to heating of the crystal lattice. In this case, if the constant 
temperature is preserved, there must be a unidirectional flux of heat from the Casimir plates to 
the heat reservoir. The existence of such an interaction between a system and a heat reservoir 
is strictly prohibited in a state of thermal equilibrium [72] and is in contradiction with its 
definition [73]. According to this definition, in thermal equilibrium all irreversible processes 
connected with the dissipation of energy are terminated. Specifically, in thermal equilibrium 
any nonzero gradients of charge carrier density and any diffusion are impossible. Thus, the 
inclusion of the screening effects and diffusion currents into the Lifshitz theory is in violation of 
its applicability conditions. 

7. Conclusions and discussion 

In the above, we have discussed several problems at the interface between field-theoretical 
description of the Casimir effect and experiments on measuring the Casimir force. The 
consideration of the Casimir energies and forces in ideal metal rectangular boxes leads to the 
conclusion that even when using ideal models it is important to take into account some general 
physical requirements. Thus, it is not productive to use the free energy which leads to the 
Casimir forces of quantum nature which increase with increasing size of the box. It also seems 
thermodynamically inconsistent to claim that the Casimir force acting on a piston is a well 
defined quantity, whereas the forces acting on all other faces of the box are excluded from 
consideration. The reason is that if the free energy is defined correctly (see Sec. 2), there is 
a uniquely defined pressure on all faces of the box equal to the negative derivative of the free 
energy with respect to the box volume calculated at constant temperature. 

An important tool for the comparison of experiment with theory is the proximity force 
approximation. In Sec. 3 we have discussed some inexact formulations which can be found in 
theoretical publications on this subject. We have also discussed recent achievements in quantum- 
field-theoretical approach to the calculation of the Casimir energies in terms of functional 
determinants and scattering matrices. This scientific direction has already obtained the first 
analytical results beyond the PFA. It is of great promise for many experimentally relevant 
applications of the theory. 

In Sees. 4 and 5 we tried to add clarity to the widely discussed problems of the precision 
of experiments on the Casimir effect and the agreement between experiment and theory. It 
was stressed that the precision of some independent measurements can be much higher than 
of respective theoretical computations using the values of parameters which may not be known 
precisely enough. In such cases the agreement of experiment with theory can also be not as 
good as the precision of the measurements. 

Finally, in Sec. 6 we have analyzed a recent theoretical approach to the thermal Casimir 



force taking into account the screening effects and diffusion currents. Using quantum Fermi- 
Dirac statistics and respective Thomas-Fermi screening length, we have applied this approach to 
calculate the Casimir free energy between two metallic plates. It was shown that the obtained 
free energy results in a violation of the Nernst heat theorem for metals with perfect crystal 
lattices. Thus, the approach under consideration is inconsistent with thermodynamics. The 
calculational results for the Casimir pressure in the configuration of two Au plates were compared 
with the results of the most precise experiment performed using a micromachined oscillator. It 
was shown that the theoretical predictions following from the inclusion of the screening effects 
are rejected by the experimental data at a 99.9% confidence level. The reason for the failure 
of this approach is the inclusion of irreversible diffusion processes violating thermal equilibrium 
which is the basic applicability condition of the Lifshitz theory. 

Phenomenologically, the Lifshitz theory combined with the generalized plasma-like dielectric 
permittivity provides a description of dispersion forces between metallic test bodies which is 
in agreement with thermodynamics and consistent with all available experimental information. 
For now there is no other theoretical approach to the thermal Casimir force between metals 
which would satisfy the requirements of thermodynamics and be simultaneously consistent with 
all measurement data. 
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